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Motivated by CoNb206 (belonging to the columbite family of minerals), we theoretically study 
the physics of quantum ferromagnetic Ising chains coupled anti-ferromagnetically on a triangular 
lattice in the plane perpendicular to the chain direction. We combine exact solutions of the chain 
physics with perturbative approximations for the transverse couplings. When the triangular lattice 
has an isosceles distortion (which occurs in the real material), the T — phase diagram is rich with 
five different states of matter: ferrimagnetic, Neel, anti-ferromagnetic, paramagnetic and incommen- 
surate phases, separated by quantum phase transitions. Implications of our results to experiments 
on CoNb206 are discussed. 

Two of the most exciting themes in condensed matter physics are the destruction of magnetic order at absolute 
zero temperature, i.e. quantum criticality pQ, and quantum fluctuations in geometrically frustrated magnets. Each 
of these themes has its archetype: the quantum Ising chain in a transverse field [T] for quantum criticality, and the 
triangular lattice Ising model for geometrical frustration [2 . Remarkably both these ideals are combined in CoNb2 06, 
a material which may be understood as a collection of weakly coupled Ising chains that sit on a triangular lattice in 
the plane perpendicular to the chain direction. Very recently, CoNb206 has been the subject of beautiful single crystal 
inelastic neutron scattering experiments by Coldea et al [3]. The observed spectrum has been shown to reflect some 
very subtle properties of one-dimensional Ising quantum critical physics, that were discovered three decades after the 
study of the quantum Ising chain was initiated [4], in a triumph of theoretical physics led by Zamolodchikov [5]. A 
full interpretation of the experimental data requires an understanding of the interplay of quantum criticality in the 
chain and the physics of frustration in the plane perpendicular to the chain. In this paper, we combine the remarkable 
insights of Zamolodchikov with the novelty that frustration brings to the phase diagram of this remarkable material, 
setting up a platform to understand the experimental data in totality. Our focus here is to chart out the global phase 
diagram of CoNb206, as well as other less well studied materials that have similar structures. 

CoNb2 06 is an insulator. A strong single ion anisotropy forces the spin on the magnetically active Co 2+ (3d 7 ) 
ions to point either "up" or "down" along the local easy axis, realizing a two- level system. The Co 2+ ions interact 
with each other through Ising interactions of the form S?Sj, where Si are the usual spin- 1/2 operators acting on the 
two level system located at site i. A magnetic field applied along the b-axis is transverse to all the easy axes and 
hence may be modeled by a term, h±Sf, which acts like a transverse field on the Ising spins. In the absence of the 
transverse field the exchange interactions will cause the Ising spins to order magnetically. On the other hand when h± 
is much larger than the exchange interactions, the Ising spins will all lie in the quantum superposition, (| T) + | l))/V% 
highlighting the fact that the transverse field is a quantum mechanical parameter with no classical analogue. The 
central question of interest that we will address in this publication is: What is the nature of the magnetic ordering 
and what is the sequence of T = phases between the two extreme limiting cases as h± is tuned? 

We begin by presenting a minimal theoretical model for CoNb2 06 which incorporates details of the columbite 
structure and which we expect to reproduce most of the features of the experiments. The lattice consists of zig-zag 
chains of Co 2+ ions with a strong ferromagnetic Ising interaction, Jo < along the c-direction. The chains form an 
isosceles triangular lattice in the a-b plane, with anti-ferromagnetic coupling to neighbors of strength J\ and J 2 as 
shown in Fig. [I] 

H = Jo ^ S%. r S% +1 . r + Jl ^ ^;r^;r+ai + ^2 ^ ^;r^;r+a 2 + ^2 ^ ^;r^;r+a 3 — ^-L ^z]V> (1) 

r r r 

a i,2,3 are the three nearest neighbors on the triangular lattice illustrated in Fig. [I] Motivated by CoNb206 we 
focus on the limit when < Ji, J2 <C |Jo|, the limit of weakly coupled Ising chains. Our focus is on the effects of 
quantum fluctuations controlled by the transverse field and hence we will not address early experimental work [7] 
which studied the material in a longitudinal field. Our results with some modifications also apply to other materials 
which are known to realize Ising spin chains weakly coupled on a triangular lattice [8]. Theoretically, a study of the 
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2-dimensional triangular Ising anti-ferromagnet in a transverse field is available [9] , which corresponds to the limit of 
Jo = of our model. We are interested in the opposite limit in this paper. Were a = J2/J1 = 1, the system would 
form a perfect triangular lattice in the basal plane. It is known that the material has a < 1 and hence the coupling 
along the chains running along the x-axis is a little stronger [6]. As we shall see, the lattice anisotropy significantly 
enriches the phase diagram. 

When Ji , J2 = the model Eq. reduces to a collection of one-dimensional Ising chains in a transverse field 
(TFIC), 

HtFIC = Jo ^2 ~ ^-L ~ H ^2 ^i' ( 2 ) 

i i i 

An ft || term is included even though we study the model Eq. |TJ) with no parallel field; the reason for this will become 
clear below. The TFIC harbors a quantum critical point (between a magnetic and a paramagnetic state when fty = 
and h± = Jo/2) in the same universality class as the well studied two-dimensional Ising field theory. Indeed, the 
singular contribution to the ground state energy of the TFIC Hamiltonian, Eq. ^ can be related to the free energy 
of the Ising field theory by a simple re-scaling of variables. The scale is set by comparing correlation functions in the 
TFIC (it can be solved exactly when ftn =0, see e.g. [10]) to the Ising field theory results [11 . We find, 

w Jo f f h ± - Jq/2 c h h\\\ 

We will make use of the scaled variables h_\_ = h± J Q J /^ 2 , h\\ = and express all energies in units of Jo/2 (ch = 
V ^ el/ 4A3 / ~ ~ 0.4016, A « 1.2824 is Glaisher's constant). Expansions for the function £ift(^_l 5 ^||) are available [11] 
in different limits for both finite h± and fty, when £ = fty/ft" <C 1 and when the magnitude of rj = — h±/hl' a is small, 

4£°(h±,h\\) = lo S ~ h l + ~ h l ( G ^ 2 + G ^ 4 + G eC 6 • • • ) - 



?h±<0 
-IFT 



8tt 

(ft ± ,ft,|) = ^logfti + fti((5 1 |£|+(5 2 |£| 2 + ...), 

^Ui;2 ,1: 16/15 



£lFT(h±,h\\) = — logft ± + ft 



8tt 



l] log T] 2 + ($0 + $1?7 + ^2^ 2 



(4) 



where a = 15/8 is a critical index of the Ising universality class. The numbers G ni G n and <I> n are universal numbers 
which are known to great accuracy, a few are reproduced in Table [T] 

In order to expand around the non-trivial TFIC limit we construct a variational wave- function, = ]J i |ft^), as a 
product of TFIC ground states, where the longitudinal field ftn = hi defining the TFIC Hamiltonian for each chain is 
taken as a variational parameters. The variational energy for our model is given by, 

^mf = ^2Si FT (h±,hi) + hjffij + ^^m~fhj, (5) 

i 

where the magnetization on the i th chain, fhi = — dSiYT / dhi. Optimization over ft^, gives hi — —Jijfhj 

_ jr 1 _ 

Emf = ^£lFT(ft_L,fti) - ^hi^^-hj. (6) 

i 

Here the labels i,j span sites in a plane (or corresponding chains), and J^ is the interaction between a pair 
of spins. In addition we have defined the rescaled matrix J^- = c\ Jij/( Jo/2), whose non-zero entries are 
J\,2 = c\ Ji,2/( Jo/2). Sometimes we prefer to work with the mean-field magnetization, fhi = J^hj, in which 

case E m f = J2 i Siyt (h±, —Jijfhj) — ^fhi^ffhj. The philosophy of our approach is similar to "chain mean-field" 
theories that has been applied to a number of unfrustrated problems [12] [13] [14] . 

We now study the quantum phase diagram (phases and phase transitions) that arise by extremizing the variational 
energy Eq. Q using the exact expansions Eq. Q for the TFIC. We present our results for the phase diagram as a 
functions of J\ and hj_, for a fixed a. 

Perfect Triangles: When J2 = J\ (i.e. a = 1), we have perfect equilateral triangles in the basal plane. Treating this 
coupling in perturbation theory, we approach the magnetic ordering from each of the two phases of the Ising chain, 
the paramagnetic phase for h± > and the magnetic phase for h± < 0. 
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When h± < 0, each isolated chain has a two-fold degenerate magnetic state. Once the chains are coupled this 
degeneracy is lifted by the exchange fields from the neighbors. Working perturbatively, we expand the energy Eq. ([6| 
order by order in J\ till the degeneracy is fully lifted. At leading order after minimizing the energy with respect to 

the magnitude of m^, E m f w fhi^ffhj, (with \fhi\ fixed) i.e. we obtain the classical triangluar lattice Ising model for 
the signs of m^. It is well known that there is an extensive degeneracy of configurations which minimize the energy [2 , 
i.e. ground states consist of all configurations with exactly one pair of parallel spins (fhi of the same sign) on each 
elementary triangle. We must continue expanding in J\ to find how this entropy is finally lifted. At next order, 
the free energy receives a contribution, SE m f « G2 \h±\ 2 ~ 2a (Jijrrij) 2 . Since G2 < 0, this perturbation prefers to 
have spins on all elementary hexagons parallel. So finding the quantum ground state of our model when h± < is 
reduced to finding the classical ground state of the model, H = J^luj) a i a j ~ K^2i (Sa°*+ e a) 2 wnere e a are all 
the neighbors of a site, when K,J > 0,X « J and o~i = ±1. We find that the ground state (proven in |I A[ ) is a 
three-sublattice ferrimagnetic state with spins pointing up on the A sub-lattice and the spins pointing down on the 
B and C sub-lattices (there are a total of 6 symmetry-related states). We denote this state as FR, see Fig. [2] 

When h± > 0, the isolated chains are in a non-degenerate gapped paramagnetic (PM) state. Perturbation theory 
in the coupling J\ leaves this state stable and hence there is always a region of PM for small enough J\. We can 
now ask if we increase J\ what magnetic order will first condense. The primary condensation takes place at (47r/3, 0). 
Since the order parameters are small, in this regime we can perform a Landau expansion for the order parameter. We 
have to consider both the complex order parameter at Q c = (47r/3, 0) which we call <£> and the real uniform component 
at (0, 0) which we call m: hi = $e l( ^ c * r + $*e _l( ^ c ' r + m. Keeping only the terms relevant to our study from the most 
general expansion for these two fields, we arrive at the following Landau continuum expression for the energy density 
E m f « J cPrei,, ignoring gradient terms, 

e L = c^|$| 2 + /^|$| 4 + a m m 2 + A 3 m($ 3 + $* 3 ) + A 6 ($ 6 + $* 6 ) (7) 

The coefficients are fixed by the parameters, G2- ) G^GQ 1 h^ and J. As J\ is increased, the condensation of takes 
place when a$ = (note however a m remains non-zero) in Eq. (JtJ) , when Jf = ^^h 2 _^~ 2 . We can then optimize m 

for a given configuration of <£, m = — Xs/(2a m )(^ 3 + <£* 3 ) and eliminate it from the energy. This shifts the A6 term 

9 3G G 8G^ ~ 2 fi 

to A6 — X^/(4a m ), which evaluates to — — ^h^ a > at the critical point. This favors states in which <£ b is 

real and negative, so that the phase defined by = \<&\e l ° is pinned to it/ 2. Magnetic order with <I> = l^e* 71 "/ 2 and 
m = corresponds to an anti- ferromagnetic state, we call AF, illustrated in Fig.[2](b). This establishes that the AF 
order forms first when the critical Jf is crossed. 

The pattern of magnetic ordering close to h± = can be found using the /^-expansion. In this regime, we find 
with 77 = that the FR state has the lowest energy and this state can be smoothly connected to the state found for 
h ± < 0. 

Putting together the results from the three expansions, when h± > 0, h± < and h± — 0, we find that the simplest 
form of the phase diagram consistent with these expansions is the one shown in Fig. |3|a). On general grounds, 
the transition between FR and AF can either be direct and and first-order, or it may proceed through a region of 
"coexistence" corresponding to states with complex <£> 6 . 

Isosceles Triangles: Now we adapt and extend the above analysis to the more complicated case when a = J2/J1 7^ 1, 
resulting from an isosceles distortion of the triangular lattice in the basal plane (as occurs in the real material). 

When h± < 0, once again there is a macroscopic degeneracy associated with the spontaneous symmetry breaking 
on each chain which will be lifted by the interactions. Exactly as in the a = 1 case, at leading order, we have 

an Ising model for m^, but now on an anisotropic lattice, E m f — fhi^-fhj. For a > 1 the ground state is a two- fold 
degenerate Neel state (Nl), shown in Fig. [2] For a < 1 (the case in CoNb20e), the situation is more complicated, and 
at leading order the ground state still has a degeneracy. The triangular lattice in the basal plane can be divided into 
"basal" chains (not to be confused with the c-axis Ising chains), see Fig.JT] each of which is Neel ordered, however each 
basal chain can be shifted by a unit without lifting the degeneracy. In fact, in mean field theory this sub-extensive 
degeneracy is not lifted at all, an artifact of this approximation. Instead, a full perturbative expansion splits the 



degeneracy at fourth order (see Appendix IB), by creating a ferro-magnetic interaction between alternate basal Neel 
chains (even and odd basal chains get locked separately). This effective interaction is a quantum fluctuation effect, not 
present in the classical h± = limit. The phase shift between even and odd basal chains is determined spontaneously 
and leads to an extra 2-fold degeneracy, resulting in a 4-fold degeneracy for the N2 state. For a fixed h± < and 
1 — a <C 1, for arbitrarily small Ji, we have shown that the N2 state is selected. If J\ is made large enough however 



one expects a transition back into the FR state. The phase boundary between N2 and FR is first order, and can be 
estimated by asking how big J\ should be so that the SE m f w G2 Y^i \h±\ 2 ~ 2a (jij m j) 2 (derived above for a = 1) term 

that picks the FR state is large enough to close the gap that the anisotropy opens; we find J^ 2_FR = }t ^ . 

^|Cr2 I {yi--\-20i) o) 
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When h± > 0, we proceed as in the case of perfect triangles. For small J\ the PM state is stable. However, at 
sufficiently strong coupling there is again a condensate but it occurs at an incommensurate wave- vector. The onset 
wave-vector is determined completely by a, it is (g*,0), where <?* = 2 cos -1 (— f ), shifted from (47r/3,0). The phase 
transition will occur again when the coefficient of the quadratic term of the (gj,0) boson mode changes sign, this 

h 2a ~ 2 

determines the critical coupling: Jf = 2 \G 2 t(2+a 2 ) ' ^ e ca ^ tms P nase w ^ n incommensurate magnetic ordering IC. 

As J i is made even stronger the system will eventually prefer to lock-in to a commensurate ordering vector. For 

1 — a <C 1, this state will be the AF state discussed above. To understand the commensurate-incommensurate phase 

transition, it is convenient to begin in the commensurate Q c = (47r/3,0) AF phase. We can then understand the 

transition to the incommensurate state as the formation of a dilute condensate of domain walls (which separate 

the 6 degenerate minima). We begin by evaluating the mean- field energy, for configurations of the form hi = 
|$| [ e i(QcTi+0(ri)) +cx _j + m? where 0( r j 

is a smooth field which creates a domain wall when it shifts between 
domains, = (27r/6)n. Owing to the anisotropy of the problem, we may assume that (apart from fluctuations), 
0(r) depends only upon x. After coarse-graining and extremizing with respect to m, we then obtain the famous 
sine-Gordon model [15], 



~ AyZ I 



dx 



L2 



{djf + S x d x 6 + Acos(n sg 0) 



(8) 



where A yz is the sample area in the y — z plane, and the couplings of the SG model can be related to the parameters in 



Eq. ([6| as worked out explicitly in Appendix I C n sg = 6. From these expressions, we can calculate the phase boundary 
between the commensurate AF state and the IC state at leading order in (1 — a): Jp IT = Jf + AciTh 2 ^ -2 (1 — a), 
where Acit = ° 47r = is a pure number. 

-2G 2V /16G1-6G 2 G 6 

Our arguments above have predicted a phase diagram shown in Fig. [3] when 1 — a <C 1. As expected from scaling 
arguments, all the phase boundaries scale as J « Ah 2 ^~ 2 when the TFIC's are close to their individual critical points, 
where A is an amplitude which depends on the transition being considered. We have succeeded in computing the 
amplitudes, A in some rather non-trivial approximations which we have detailed in the body of the work. Since a < 1 
in CoNb206, our results should apply to this material. Let us now turn to a discussion of the experimental response 
of each the phases as h± is increased from zero and to the nature of the phase transitions between them. 

Neel phase: This is the phase N2 (illustrated in Fig. |2| which appears with weak external fields. It has no net 
magnetization and exhibits an elastic neutron peak at momentum (f,0) ((0, ^,0) in the orthorhombic co-ordinate 
system). This is in agreement with zero field experiments [6 , 7 , though further-neighbor interactions not present in 
our model could play a role at h± =0. 

Ferri-magnetic phase: As h± is increased, we expect a first order phase transition into the FR state. The transition 
N2-FR is driven by the increase of quantum fluctuations of the Ising spins. The FR phase has a net uniform magnetic 
moment and thus is characterized by elastic neutron peaks both at zero momentum and at the commensurate ordering 
vector (47r/3,0). As h± is increased we expect the FR order to weaken causing a reduction of the Bragg peaks; the 
uniform magnetization m and hence zero momentum peak goes to zero as h± is increased much faster (like the third 
power) than the peak at (47r/3,0). 

Anti- Ferromagnetic phase: The AF state also orders at (47r/3,0), but it has no uniform component of the magne- 
tization in the z-direction and is hence analogous to an anti-ferromagnet. It is entirely different from the N2 state 
however, and is truly a result of the geometric frustration present. The difference between the AF and FR states 
could be probed by tilting the field in the b — c plane - one expects a weak first order transition at zero tilt in the 
FR case but not the AF one. 

Incommensurate: As h± is increased further still, there is a transition into the IC phase that is driven by a 
condensation of domain walls in the AF phase. The IC phase is characterized by a dense set of Bragg peaks and a 
massless goldstone (phason) mode. The primary ordering vector, q x , shifts continuously in this phase from q x = 4tt/3 
in t he AF phase to value q* = 2cos _1 (— a/2) at the transition to the PM, the evolution of q x (h±) is described in 

ED- 

Paramagnetic: Finally when the quantum fluctuations due to h± are strong enough, magnetic ordering is completely 
suppressed and the system enters a quantum paramagnetic phase. This state is gapped and breaks no symmetries of 
the model, Eq. 0. 

Quantum Phase Transitions: There are four quantum phase transitions that we have predicted in the evolution of 
phases as h± is tuned. Although the location of these transitions relies on mean field theory, their critical properties 
are universal. The first transition between N2-FR is expected to be strongly first-order. The second transition is 
between the FR-AF states. It can take place either by a direct first order transition or through a coexistence of the 
two phases. Within mean field theory, we find the latter, though this conclusion may be sensitive to the approximation, 
and to small perturbations of the Hamiltonian. The third transition between the AF-IC state is continuous and is 
of the commensurate-incommensurate type. We have calculated the leading singularity of the ordering vector here. 
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The final and fourth transition is between IC-PM, which we expect to be in the XY universality class in 4 classical 
dimensions. Since this is at the upper critical dimension we expect mean field exponents with log corrections. The 
transverse magnetization as a function of the external field, h±, should provide a good method for locating the phase 
transitions. The magnetization itself should have discontinuities at the first order transitions and it's first derivative 
should have singularities at the locations of the second order transition. 

In summary, we have studied the model, Eq. 0, for CoNb206 in an external magnetic field that is aligned along the 
crystallographic b-axis, which results in a fully quantum mechanical "transverse field" . Application of the transverse 
field drives the frustrated Ising model under study into a number of fascinating phases through quantum phase 
transitions. We have outlined the basic properties of each of these phases here. Detailed calculations of the excitation 
spectrum, measurable by neutron scattering, are also possible, and will be presented in a follow-up publication. An 
interesting question we leave for future exploration is the nature of the finite-temperature phase diagram of this model. 

Acknowledgements.- We thank Radu Coldea for stimulating discussions and acknowledge support from the Packard 
Foundation and National Science Foundation through grants DMR-0804564 and PHY05-51164. 



I. SUPPLEMENTARY MATERIALS 
A. Ground state for a — 1 and h± < 

Here we prove that the ground state of the model on the perfect triangular lattice, H — JYluj) (7 i (7 3 ~ 

K^ /i (J2 a a i+ea) 2 i f° r J,K > and K/J <C 1 is the FR state. The proof is as follows: We need to find which 
Ising configuration has the lowest energy with the K term. Lets imagine an Ising configuration on an infinite lattice 
with /g, / 5 , /4, fs the fraction of sites that have 6, 5, 4, 3 anti-parallel neighbors. Then /6 + /5 + /4 + /3 = 1- Also a third 
of the bonds have to be frustrated (host parallel spins) in any Ising configuration which implies, f$ + 2/4 + 3/3 = 2. 
The contribution of the K term to the energy in this language is: E = —K(36fe + I6/5 + 4/4). Since we have two 
extra equations, and four unknowns, we can rewrite the equation for E completely in terms of fs and /g, obtaining 
E = — if (4+8(/3 + /e)). Since fs-\-fe < 1, the lowest energy is clearly obtained only with 3, 6 sites. This configuration 
has /e = 1/3 and fs = 2/3. The only such configuration is the ferri-magnetic state described above. 



B. Ground state for a < 1 and h± < 



In the section we show that the ground state when a < 1 and when the chains are magnetically ordered is the N2 
state for arbitrarily small inter-chain coupling. The goal here is to develop a systematic perturbative expansion about 
the limit of decoupled Ising chains. On the magnetic side there is a huge degeneracy in this limit because each chain 
can have its magnetization mi point up or down. Working in imaginary time, we can set it up as follows, 

Z f e - Sld ~ Sc , o 

— £ — 



Z fe~ 



<e-*«), (9) 



where S c = J drS^(r)^S](r) and all averages are taken over the decoupled limit. Note that in the zero temperature 
limit, the partition function is simply given by, Z = e~^ E ^ s . We can calculate the average perturbatively using the 
cumulant expansion, 

(e~ S °) = e _ [^ c ) _ A"^c)con+|[<5'c)con-^! <5'c)con + --- ] ^ {W) 

where the subscript "con" indicates that we should only keep connected averages. When we have to evaluate the 
connected averages we can factor the average into products of averages each of which acts within a single chain, since 
the chains are decoupled at th order. For an odd number of spins, these single chain averages will be replaced by, 
e.g. (Sf(r)) — > fhi\(S z (r))\ with fhi = ±1. Then we can ask which configuration of fhi produces the lowest energy 
and keep evaluating higher order contributions till any "accidental" degeneracy is completely lifted. The coeffecients 
of the interactions between fhi are determined by integrals over correlation functions of the d=l Ising chain. 

We now turn to an evaluation of this approach for the case of interest, when a < 1. At leading order, we have the 
term rhi-^-rhj, which gives Neel configurations along the stronger bonds that form a chain in the basal plane (not 
to be confused with the Ising chains in the transverse direction!), but each basal chain can be translated by a unit 
leading to a 2 7Vchain degeneracy (sub-extensive entropy). We expect that the degeneracy is lifted by locking-in of spins 
on alternate basal chains, but we still need to determine whether the locking-in is ferromagentic or anti-ferromagnetic. 
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At second order, we generate ferro two-spin interactions between alternate basal chains ex —JijJikfhjm k but it is easy 
to see that these cancel when summed over the sites of the basal chain due to the Neel order, i.e., the mean field from 
all the sites in a chain add up to zero. The degeneracy is not lifted at this order. At third order, we generate both 
two-spin and four-spin couplings. It turns out that both the two- and four-spin couplings cancel out. At fourth order, 
we have a large number of generated two-spin couplings and four-spin couplings. Turns out that the four-spin coupling 
does not lift the degeneracy, but a careful summation of the 86 two-spin couplings does! The two-spin interactions at 
fourth order are ferromagnetic and hence alternate chains are locked ferromagnetically as shown in Fig. [2] 



C. Commensurate-Incommensurate Transition 



In this section we calculate the energy of a domain wall in the commensurate order parameter as a function of the 
external parameters. The simplest commensurate-incommensurate transition (CIT) occurs when single domain walls 
condense [15], i.e. when the energy of a single domain wall goes to zero. We can estimate this by a calculation of the 
domain wall energy. A domain wall corresponds to a spatial profile of the phase 0(r) of the complex order parameter 
<I> introduced earlier in the context of Eq. We hence consider the following form for the field hi at site ri, 



e *(QcT<+*(r,)) +c . c . 



m. (11) 



Here, Q c is a commensurate wave vector (47r/3,0) and 0(r) is a slowly varying functions of r. We substitute this 
expression into our variational energy, Eq. (pi and then derive a coarse grained energy density for the field 0(r). Next 
we integrate out m, like we discussed below Eq. ([7]), arriving at the Sine-Gordon long- wavelength energy for 0, Eq. (J8|, 
where we have: 

_ -8^4 + ^2^672-60 

A 6 - ^ ">± j 

A3 = —SG^h^ 4a , 
a m = -3G 2 h 2 - 2a , 

A = 2W 6 ( X f A ' 



= -\$\ 2 G 2 h 2 ± 



4a n 

-4a 



5 X = -m 2 2G ^- a \ i- 2 *. 

Except for 5 X , we evaluated all these couplings with a = 1. 

Now we are all set to calculate the transition from the state with commensurate (47r/3,0) ordering to the incom- 
mensurate phase. Since the 5 term has only an x gradient, we will look for solutions that are translationally invariant 
in the y— direction, dropping all derivatives in this direction. Noting that S only gives a boundary term, it may 
expressed simply as the winding number N w = [6(x — L) — 6(x = 0)] , the number of times the field moves from 
one its six minima to an adjacent one. The solution for N w = ±1 is well known, 

0(x) = - tmr\e ±n ^ {x - Xo) ), (12) 
n 

where n sg = 6 in our case. 

The energy of a single domain wall is positive for \8\ < \5 C \ = 4y ' K x \ji\, and the wave vector is commensurate. For 
1 5 1 > \S C \, domain walls condense and the wave vector is incommensurate. In order to compute the phase boundary, 
we need the magnetization for a given h± and J\ close to the transition when a = 1. Within our chain-MFT we can 
compute this by optimizing the mean-field energy. We find, 

r 4a-2 j _ jc 

where Jf = h 2 *~ 2 /{<o\G 2 \). Now using this as input we can calculate for 1 — a <C 1 the commensurate-incommensurate 
phase boundary. We find at leading order in (1 — a): 

J? lT = Jl + AciT~h 2 r\ l - ") ( 14 ) 
where Aqyv = G47r = is a pure number. 
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FIG. 1: Schematic structure of Co 2+ ions in C0ND2O6. Ising spins localized on each Co 2+ form ferromagnetic chains running 
along the c-direction; Jo is the strength of the exchange along these c-axis chains. In the basal a-b plane perpendicular to the 
c-axis, the chains form a triangular lattice. The stronger anti-ferro J\ bonds parallel to the b-axis divide the basal triangular 
lattice into "basal chains" (marked with thick black lines), the weaker anti-ferro J2 = aJi bonds (marked with thin green lines) 
connect neighboring "basal chains" . 1 — a <C 1 in CoNb206 is a measure of the departure from a perfect basal triangular lattice 
in this material. 



D. Ordering vector in the IC 

Once the hi condense, initially at some incommensurate vector <?*, it is clear that the ordering vector will change 
continuously in the IC phase as h± is tuned. We can calculate the evolution of the ordering vector from q* by 
performing a harmonic expansion valid close to the transition. Approaching the IC phase at high h± from the PM 
state, we look for mean-field solutions including a higher harmonic, hi = ^^e lk ' ri + ^^e l ^ k ' ri + ex., which minimize 
the energy, Eq. We can write the energy density as, e = otk\&k\ 2 + Uk\&k\ 4 + Oisk\^3k\ 2 + \k(®k®3k + c.c). 

a 2 A 2 / \ ^ 

Optimizing over an d $3/e, we find that minimum energy density is e = — ^ — ( ) . If = then the 

optimal k — q*, but in practice there is a shift, that to leading order in (1 — a) and h± — h\ is given as, 

32(2o-2)% f h ± -h c ± \ 2 

S Qx = Tlx ( X - a ) ( J ( 15 ) 



9^ V h' 



from which we conclude that the shift is positive for (1 — a) <C 1. 

We now turn to the opposite edge of the IC phase, when it borders a commensurate phase, i.e. the AF phase 
where the ordering takes place at (47r/3, 0). Once we cross into the IC phase the finite density of domain walls shifts 
the ordering vector. The shift is given by 5q x = 27rn^, where is the linear density of domain walls. In order to 
estimate n^, we write an expression for the energy of a lattice of domain walls including an interaction energy between 
domain wall pairs from the SG model (assuming n^w <C 1), 



E sg (n d ) ( SVkA 2ttS x 



A yz 



n d + Un d e wn * , (16) 



where U > is the strength of repulsion between domain walls and w = ^—^/j is the width of domain walls and 

controls the range over which they interact. When S x > S c = 4V^A/tt, there is a finite density of domain walls, 
n d ^ w\n\5 -s I • F rom this we can estimate the leading singularity in the shift in the ordering vector close to the CIT, 

-27T 

Sq * = ~1 — it ^cItT ( 17 ) 

w log[h± — hj_ L \ 



where w is the width of the domain wall at the CIT. From our derived expressions Eq. (12) for the coupling constants 
in the SG model, we can calculate the leading dependence of the domain wall width at the CIT, we find, w = ^(i-^n ' 
A plot of the two asymptotes for the ordering vector is shown in Fig. [4] 
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(a)Ferri-magnetic (FR) 



(b)Anti-Ferromagnet (AF) 




(c)Neel 1 (Nl) 



MM 
M MM 
MM r 1 
MM M 

(d)Neel 2 (N2) 

FIG. 2: Magnetic ordering patterns of FR, AF, Nl and N2 phases, (a) Dividing the lattice into A, B and C sub-lattices 
(marked red, blue and black), the FR phase corresponds to putting the spins on two of these parallel and the spins on the 
third anti-parallel, leading to a six- fold degeneracy, (b) The AF state is also a three sub-lattice state but with the spin (Sf) 
on one site fluctuating to zero (c) the two-fold degenerate Nl phase, obtained when a > 1, can be thought of most simply as 
the bi-partite (red and blue) Neel phase obtained when the Ji bonds are set to zero, (d) the four- fold degenerate N2 phase 
obtained when a < 1 corresponds to a Neel ordering along the basal chains ( Ji bonds) and a ferromagnetic locking of alternate 
chains (shown by dashed lines). 





n = 1 


n = 2 


n — 3 


G'2n 


-1.84522 


8.33370 


-95.1689 


G n 


-1.35783 


-0.0489532 


0.0387529 


$n-l 


-1.19773 


-0.318819 


0.1108915 



TABLE I: Leading coefficients of the energy expansion of the Ising field theory, appearing in Eq. Q, following Ref. 1111 
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FIG. 4: Primary ordering vector in the IC phase. We have computed (in Section ID) the asymptotics expanding around 
the CIT and the transition from the PM phase. The evolution of the q x is very sharp as shown. The inset shows the same 
data on a linear- log plot. The plot has been made for a = 0.99 and Ji =0.1 which gives, q% w 4.17726, h° ± w 1.05586 and 
/i? IT w 0.99842. 
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